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Abstract 



We overview numerous algorithms in computational D-module theory together 
with the theoretical background as well as the implementation in the computer 
algebra system Singular. We discuss new approaches to the computation of Bern- 
stein operators, of logarithmic annihilator of a polynomial, of annihilators of rational 
functions as well as complex powers of polynomials. We analyze algorithms for local 
Bernstein-Sato polynomials and also algorithms, recovering any kind of Bernstein- 
Sato polynomial from partial knowledge of its roots. We address a novel way to 
compute the Bernstein-Sato polynomial for an affine variety algorithmically. All 
the carefully selected nontrivial examples, which we present, have been computed 
with our implementation. We address such applications as the computation of a 
zeta-function for certain integrals and revealing the algebraic dependence between 
If*. | pairwise commuting elements. 
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1 Introduction 

Constructive D-module theory has been dynamically developing throughout the last years. 
There are new approaches, algorithms, implementations and applications. Our work 
on the implementation of procedures for D-modules started in 2003, motivated among 
other factors by challenging elimination problems in non-commutative algebras, which 
appear e. g. in algorithms for computing Bernstein-Sato polynomials. We reported on 
solving several challenges in [20]. A non-commutative subsystem Singular: Plural [14] 
of the computer algebra system Singular provides a user with possibilities to compute 
numerous Grobner bases-based procedures in a wide class of non-commutative G-algebras 
[22]. It was natural to use this functionality in the context of computational .D-module 
theory. Nowadays we present a D-module suite in Singular consisting of the libraries 
dmod . lib, dmodapp . lib, dmodvar . lib and bf un . lib. There are many useful and flexible 
procedures for various aspects of D-module theory. These libraries are freely distributed 
together with Singular [11]. 

There are several implementations of algorithms for D-modules, namely the experi- 
mental program kan/smI by N. Takayama [38], the bf ct package in Risa/Asir [31] by 



M. Noro [30] and the package Dmodules .m2 in MACAULAY2 by A. Leykin and H. Tsai [40]. 
We aim at creating a -D-module suite, which will combine flexibility and rich functionality 
with high performance, being able to treat more complicated examples. 

In this paper we do not present any comparison between different computer algebra 
systems in the realm of -D-modules, referring to [20] and [1]. However, comparison in the 
latter articles shows, that our implementation is superior to kan/smI and Macaulay2 
and in many cases more powerful than Risa/Asir. 

Here is the list of problems we address in this paper: 

• s-parametric annihilator of / (Section 3, see also [20, 1]), 

• annihilator of f a for a G C (Section 4, see also [34]), 

• annihilator of a polynomial function / and of a rational function f/g (Section 4), 

• 5-function with respect to weights for an ideal (Section 5, see also [1]), 

• global and local Bernstein-Sato polynomials of / (Section 6), 

• partial knowledge of Bernstein-Sato polynomial (Section 6.4, see also [20]), 

• Bernstein operator of / (Section 7), 

• logarithmic annihilator of / (Section 8), 

• Bernstein-Sato ideals for / = f\ ■ . . . ■ f m (Section 9, see also [20]), 

• annihilator and Bernstein-Sato polynomial for a variety (Section 10, see also [1]). 

We describe both theoretical and implementational aspects of the problems above and 
illustrate them with carefully selected nontrivial examples, computed with our implemen- 
tation in Singular. In Section 3.3, we give yet another alternative proof for the algorithm 
by Briangon-Maisonobe for computing AnnD n r s i(/ s ), presented in [1]. Notably, this de- 
livers additional structural information. In Section 7, we compare several approaches for 
the computation of Bernstein operators. Using the method of principal intersection, we 
formalize several methods for computing Bernstein-Sato polynomials. Following Budur 
et al. [7] and [1], we report on the implementation of two methods for the computa- 
tion of Bernstein-Sato polynomials for affine varieties in a framework, which is a natural 
generalization of our approach to the algorithm by Briangon-Maisonobe. 

Notations. Throughout the article IK is assumed to be a field of characteristic zero. By 
R we denote the polynomial ring K[xi, . . . , x n ] and by / G R a non-constant polynomial. 

We consider the n-th Weyl algebra as the algebra of linear partial differential operators 
with polynomials coefficients. That is D n = D(R) = K(xi, . . . , x n , d\, . . . , d n \ {diXi = 
Xidi + 1, diXj = Xjdi, i 7^ j}). We denote by D n [s] = D{R) ®k K[si, . . . ,s n ] and drop the 
index n depending on the context. 

The ring R is a natural -D n (i?)-module with the action 



Xi • f(xi, ...,x n )=Xi- f(x u ...,x n ), di* f(x u 



df(x 1 ,...,x n ) 
dx. 



Working with monomial orderings in elimination, we use the notation x ^> y for u x is 
greater than any power of y" . 

Given an associative IK-algebra A and some monomial well-ordering on A, we denote 
by lm(/) (resp. lc(/)) the leading monomial (resp. the leading coefficient) of / G A. 
Given a left Grobner basis G C A and / G A, we denote by NF(/, G) the normal form 
of / with respect to the left ideal a(G). We also use the shorthand notation h — >h f 
(and h —¥ f, if H is clear from the context) for the reduction of h G A to / G A with 
respect to the set H. If not specified, under ideal we mean left ideal. For a,b G A, we 
use the Lie bracket notation [a,b] := ab — ba as well as the skew Lie bracket notation 
[a, b] k := ab - k ■ ba for k G K*. 

It is convenient to treat the algebras we deal with in a bigger framework of G-algebras 
of Lie type. 

Definition 1.1. Let A be the quotient of the free associative algebra K(xi, . . . ,x n ) by 
the two-sided ideal J, generated by the finite set {xjXi — XiXj — dij} VI < i < j < n, 
where d^ G K[x\, . . . , x n ]. A is called a G-algebra of Lie type [22], if 

• \/l<i<j<k<n the expression dijXk — Xkdij + Xjdik — dikXj + djkXi — Xidjk reduces 
to zero modulo / and, 

• there exists a monomial ordering -< on K.[x\, . . . , x n ], such that lm(dy) -< XiXj, Vi < j. 

G-algebras are also known as algebras of solvable type [17, 23] and PBW algebras [8]. 
We often use the following. 

Lemma 1.2 (Generalized Product Criterion, [22]). Let A be a G-algebra of Lie type and 
f,g G A. Suppose lm(f) andlm(g) have no common factors, then spoly(/, g) — >{f, g \ [fi9\- 

2 Challenges for Grobner bases engines of Singular 

Since the very beginning of implementation of algorithms for D-modules in SINGULAR 
there have been intensive interaction with the developers of Singular. Numerous chal- 
lenging examples and open problems from constructive D-module theory were approached 
both on the level of libraries and in the kernel of Singular and Singular: Plural. 
This resulted in several enhancements in kernel procedures and, among other, motivated 
M. Brickenstein to develop and implement the generalization of his slimgb [6] (slim 
Grobner basis) algorithm to non-commutative G-algebras. Indeed slimgb is a variant 
of Buchberger's algorithm. It is designed to keep polynomials slim, that is short with 
small coefficients. The algorithm features parallel reductions and a strategy to minimize 
the weighted lengths of polynomials. A weighted length function of a polynomial can be 
seen as measure for the intermediate expression swell and it can consider not only the 
number of terms in a polynomial, but also their coefficients and degrees. Considering the 
degrees of the terms inside the polynomials, slimgb can often directly (that is, without 
using Grobner Walk or similar algorithms) compute Grobner bases with respect to e. g. 
elimination orderings. The procedure slimgb demonstrated very good performance on 
examples from the realm of D-modules [20] , which require computations with elimination 
orderings. 



As it will be seen in the paper, various computational questions, arising in D-module 
theory, use much more than Grobner bases only. Among other, a transformation matrix 
between two bases (called Lift in [15]), the kernel of a module homomorphism (called 
Modulo in [15]) and so on must be applied for complicated examples. On the other 
hand, the standard std routine for Grobner bases, generalized to non- commutative G- 
algebras, is used together with slimgb for a variety of problems. Since the beginning of 
development of the .D-module suite in Singular, these functions have been enhanced: 
they became much faster and more flexible. The effect of the use of the generalized Chain 
Criterion (cf. [15]) in Grobner engines is even bigger in the non-commutative case, due 
to the discarding of multiplications, which complexity is increased, compared with the 
commutative case. On the contrary, the generalized Product Criterion (Lemma 1.2) plays 
a minor role in the implementation, since the complete discarding of a pair generalizes to 
the computation of a Lie bracket of the pair members. 

The concept of ring list, introduced in Singular in 2004, enormously simplified the 
process of creation and modification of rings (like changing the monomial module ordering, 
regrouping of variables, modifying non-commutative relations, working with parameters 
of the ground field etc.). Especially in the D-module setting we modify rings often, create 
a new one from existing rings and equip a new ring with a new ordering. Thus, with ring 
lists the development of such procedures became much easier and the corresponding code 
became much more manageable. 

We have to mention, that in the meantime the implementation of non-commutative 
multiplication in the kernel of Singular: Plural has been improved as well. 

3 s- parametric annihilator of / 

Recall Malgrange's construction for / = f\ • ■ • f p G K[xi, . . . ,x n ]. Consider the algebra 
W p+n , being the (p + n)-th Weyl algebra 

K({tj,dtj \l<j<p}\ {[dtj,t k ] = S jk }} (g) K K({xi,di | 1 < i < n} | {[d h x k ] = 5 ik }}. 

Moreover, consider the left ideal in W p+n , called Malgrange ideal 

p at 

If'=({ tj ~ fj, E 7T dt J +di I !< J < P» !<*<«} )■ 

j=1 OXi 

Then for s — (si, . . . , s p ) we denote f s := fl 1 ■ ■ ■ fp p . Let us compute 

IfDKlitjdtj^ixi^Xi | [di,Xi] = 1) C D n [{tjdtj}) 

and furthermore, replace tjdtj with — Sj — 1. The result is known (e. g. [34]) to coincide 

with the s-parametric annihilator of f s , Ann Dn ^(f s ) = {Q(x, d, s) G D n [s] \ Q • f s = 0}. 

There exist several methods for the computation of s-parametric annihilator of f s . 



3.1 Oaku and Takayama 

The algorithm by Oaku and Takayama [32, 34] was developed in a wider context and 
uses homogenization. Consider the K-algebras T := K[ti, . . . ,t p ], D' := D(T) and H : = 
D n cg) K D' p ® K K[ui, . . . , u p , Vi, . . . , v p ]. Moreover, let / below be the (u, v )-homogenized 
Malgrange ideal, that is the left ideal in H 



i> 



df, 



I = ( {tj - Ujfj, J2 -f£ U * dt 3 + ^' U 3 V 3 ~ 1 })- 



k=l 



Oaku and Takayama proved, that Anno n [ s ]{f s ) can be obtained in two steps. At first 
{uj,Vj} are eliminated from / with the help of Grobner bases, thus yielding I' — I D 
(D n ® K D'). Then, one calculates /' fl (D n ® K K[{— tjdtj — 1}]) and substitutes every 
appearance of tjdtj by — Sj — 1 in the latter. The result is then Ann Dn ^(f s ). 

3.2 Briangon and Maisonobe 

Consider S p = ~K{{dtj,Sj} \ dtjSk = s^dtj — Sj^dtj) (the p-th shift algebra) and S' = 
D n ®k Sp. Moreover, consider the following left ideal in S': 



i = ({8 j +f j m j ,J2^*+ d *n- 



Briangon and Maisonobe proved [5] that Ann£> n [ s ](/ S ) —In D n [s\, . . . , s p ] and hence the 
latter can be computed via the left Grobner basis with respect to an elimination ordering 
for {dtj}. 

3.3 Another alternative proof of Briancpn-Maisonobe's method 

Here we give yet another [1] computer algebraic proof for the method by Briangon and 
Maisonobe. 

Throughout this section, we assume 1 < % < n and 1 < j < p. Define 

E := K({tj, dtj,Xi, di, Sj} | {[di, Xi] = 1, [dtj,tj] = 1, [t k , Sj] = S jk tj, [dt k , Sj] = -S jk dtj}). 

Let B = D n [s] be a subalgebra of E, generated by {xi,di,Sj}. Then the Briangon- 
Maisonobe method requires to prove [1], that 

i=i * 
Theorem 3.1. Let us define the following polynomials and sets: 

P Qf 

9i := di + J2 ir dt k, G = { 9i }, T = {tj - fj}, S = { Sj + fjdtj}. 

k=i X% 

Let A be a (possibly empty) subset of {l,...,p}. Define Ma := G U {t k — f k \ k G 
A}u{ Sj + fjdtj\je{l,...,p}\A}. 



(a) For any A, the elements of Ma commute pairwise. In particular, so do G U T 
and G U S. 

(b) Consider an ordering -<, satisfying {tj} ^> {x.i\, {di,Sj} ^> {xi,dtj}. Then any 
subset of GUTU S is a left Grobner basis with respect to -< . In particular, so is the 
set Ma for any A. 

(c) The elements of Ma. are algebraically independent. 

(d) For any A, the Krull (and hence the Gel'fand-Kirillov) dimension ofK[M\] is n+p. 

(e) For any A ; K[Ma] is a maximal commutative subalgebra of E. 
Proof, (a) Computing commutators between elements, we obtain 

\9i, 9k] = dt 3 J>, |&] + dt 3 £[|&, d k ] = dtj £([*, ^] - fa, §£]) = 0, 

df- 

[tk - fk, 9i] = Yl &rf tfe ' dtj ^ ~~ ^ fc ' ^ = °' fc " f*> tk " ^ = °' 
j 

[st + fidti, Sj + fjdtj] = fj[si, dtj] - fi[sj, dti] = 0, [sj + fjdtj, g t ] = [sj, d t ]+ 

+dt j [f j , d t ] + j^ ^-[ Sj , dt k ] + [f 3 dt 3 ,J2 §^ = d A dt > ~ ft> fi\^i = °- 

fc=i * fc=i % i 

The only nonzero commutator arises from 

[tk - fk, Sj + fjdt 3 ] = [t k , Sj] + fj[t k , dtj] - [f k , Sj] - [f k: fjdtj] = 6 jk (t k - fk)- 

However, according to the definition, only one of these elements belongs to M\ for any A. 

(b) We run Buchberger's algorithm by hands. Due to the ordering property, for each 
pair the generalized Product Criterion is applicable. Hence using (a) we see, that most 
s-polynomials reduce to commutators, which are zero except spoly(t fe — f k , s 3 + fjdtj) = 
8jk(t k — f k ), which reduces to zero modulo the first polynomial. Thus, any subset including 
M\ is indeed a Grobner basis. 

(c) Using pairwise commutativity, we employ the Commutative Preimage Theorem 
from [19]. It states, that the ideal of algebraic dependencies between pairwise commuting 
elements {h k | 1 < k < m} C E can be computed as 

E Ok K[ci, . . . , Cm] D {{hi - Ci}) n K[ci, . . . , c m ], 

where q are new commutative variables, adjoint to E. In this elimination problem one 
requires an ordering on E %K[c], preferring variables of E to q's. For such an ordering, 
one needs to compute a Grobner basis. Now, take {hi} := Ma, 1 < i < p + n, and run 
Buchberger's algorithm with respect to the same ordering as in (b). Thus we are again 
in the situation, where the Product Criterion applies, hence [hi — Ci, h k — c k ] =0 since 
[hi, h k ] = by (b) and q are central. Hence, {hi — q} is a left Grobner basis and by 



the elimination property ({hi — q}) fl K[ci, . . . , c m ] = 0, that is {hi} are algebraically 
independent. 

(d) By (c), M\ generates a commutative ring with no algebraic dependence between 
its elements, so the Krull dimension is the cardinality of Ma, that is n + p. Since Ma is 
isomorphic to a commutative polynomial ring by (c), its Gel'fand-Kirillov dimension over 
the field K is n + p as well. 

(e) With respect to the ordering from (b), the leading monomials of the generators 
are {d±, . . . ,d n } U {tk \ k G A} U {sj \ j £ A}. Assume, that there exists an element 
in E \ K[Ma], which commutes with all elements in M\. Then its leading monomial 
must belong to the subalgebra F, generated by {x±, . . . ,x n } U {sk \ k G A} U {tj \ j £ 
A}u{dt±, . . . , dtp}. Since the center of E is IK, we consider centralizers of elements. Taking 
F' = r\ke\C(tk — fk) H F, we see that an element from it can have no {dt^, Sk \ k G A}. 
Considering F" = rijC(yj) fl F', we exclude {x±, . . . ,x n }. Thus we are left with the 
subalgebra, generated by F = {dtj, Sj | j G: A}. But no element of it can commute with 
{sj + fjdtj | j $. A} except constants. Hence the claim. □ 

We want to eliminate both {tj} and {dtj} from an ideal, generated by G U S U T. 
By using an elimination ordering for {tj} we proved in (b) above, that G U S U T is a 
Grobner basis. Hence, the elimination ideal is generated by G U S and we can proceed 
with eliminating {dtj} from GUS, which is exactly the statement of Briangon-Maisonobe 
in Section 3.2. 

3.4 Implementation 

We use the following acronyms in addressing functions in the implementation: OT for 
Oaku and Takayama, LOT for Levandovskyy's modification of Oaku and Takayama [20] 
and BM for Briangon-Maisonobe. Moreover, it is possible to specify the desired Grobner 
basis engine (std or slimgb) via an optional argument. 

For the classical situation / = /i, the procedure Sannf s(/) computing Ann Dn [ s ](/ s ) C 
D n [s] uses a "minimal user knowledge" principle and chooses one of three mentioned algo- 
rithms. Alternatively, one can call the corresponding procedures Sannf sOT, Sannf sLOT, 
Sannf sBM directly. 

For the annihilator oi f — f± • • • f p , see Section 9. 

Example 3.2. We demonstrate, how to compute the s-parametric annihilator with 
Sannf s. This procedure takes a polynomial in a commutative ring as its argument and 
returns back a Weyl algebra of the type ring containing an object of the type ideal 
called LD. Note, that the latter ideal is a set of generators and not a Grobner basis in 
general. 

LIB "dmod.lib"; 

ring r = 0,(x,y),dp; // set up the commutative ring 

poly f = x~3 + y~2 + x*y~2; // define the polynomial 

def D = Sannfs(f); setring D; // call Sannf s and change to ring D 

LD = groebner (LD) ; LD; // compute and print Groebner basis 

==> LD [1] =2*x*y*Dx-3*x~2*Dy-y~2*Dy+2*y*Dx 

==> LD [2] =2*x~2*Dx+2*x*y*Dy+2*x*Dx+3*y*Dy-6*x*s-6*s 



==> LD [3] =x~2*y*Dy+y~3*Dy-2*x~2*Dx-3*x*y*Dy-2*y~2*s+6*x*s 
==> LD [4] =x~3*Dy+x*y~2*Dy+y~2*Dy-2*x*y*s-2*y*s 

==> LD [5] =2*y~3*Dx*Dy+3*x~3*Dy~2+x*y~2*Dy~2-4*x~2*Dx~2-8*x*y*Dx*Dy-2*x~2*Dx 
-4*y~2*Dx*s+6*x*y*Dy+12*x*Dx*s-10*x*Dx-6*y*Dy+12*s 

4 Annihilators of polynomial and rational functions 

4.1 Annihilator of f a for a G C 

It is known (e. g. [34]) that for any a G C, D n / AmiD n (f a ) is a holonomic .D-module. 
In the procedure annf special from dmod.lib we follow Algorithm 5.3.15 in [34]. Given 
/ and a, we compute Anno„[ s ]{f s ) C -D n [s], the Bernstein-Sato polynomial of / (cf. 
Section 6.1) and its minimal integer root sq. Then, if a — (sq + 1) G N, according 
to Algorithm 5.3.15 in [34] we have to compute a certain syzygy module in advance. 
Otherwise, ArniD n (/ a ) = Annc n [ s ](f s ) \ s =a is obtained via substitution. 

Example 4.1. In this example we show, how one computes the annihilator of 2xy. 

LIB "dmod.lib"; option(redSB) ; option(redTail) ; 

ring r = 0,(x,y),dp; poly g = 2*x*y; 

def A = Sannfs(g); setring A; // compute Ann(g's) 

LD = groebner(LD) ; LD; // GB of the ideal Ann(g~s) 

==> LD[l]=y*Dy-s 

==> LD[2]=x*Dx-s 

def B = annf sO(LD,2*x*y) ; setring B; // compute BS polynomial 

BS; // the list of roots and multiplicities of BS polynomial 

==> [1]: 

==> _[1]=-1 

==> [2]: 

==> 2 

// so, the minimal integer root is -1 

setring A; // need to work with Ann(g's) again 

ideal I = annf special (LD,2*x*y,-l, 1) ; 

// the last argument 1 indicates that we want to compute f~l 
print (matrix(I) ) ; // condensed presentation 

==> Dy~2, y*Dy-l, Dx~2, x*Dx-l 

4.2 Alternative for an annihilator of f m 

Computing a syzygy module in the previous algorithm can be expensive. Therefore we 
note, that for a = m G N we better use an easier approach. 

Lemma 4.2. Let g G K.[xi, . . . ,x n ]. Consider the homomorphism of left D n -modules 
ip : D n -> D n /(d 1 , . . .,d n ), ip{\) = g. Then Ann Dn (g) = ker^. 

Proof. Note, that K[x±, . . . , x n ] = D n /(di, . . . , d n ) as left _D n -modules. Hence we can view 
g as the image of 1 under if). Then AniiD n (g) = {a G D n \ a • g = 0} = {a G D n \ ag G 
(<9i,...,<9 n )} = kerV>. □ 



Remark 4.3. Hence, given any element / G K[xi, . . . ,x n ], Anno n (f) can be computed 
via the kernel of a module homomorphism (algorithm Modulo) which amounts to just one 
Grobner basis computation. Moreover, it does not use elimination and hence is clearly 
more efficient in the special case g = f n for / G R, n G N, than the more general method 
in Section 4.1. Notably this method can be generalized to various other operator algebras, 
see [35] for details. The corresponding procedure in dmodapp.lib is called annPoly. 

Remark 4.4. Yet another improvement can be achieved in the computation of the min- 
imal integer root of the Bernstein-Sato polynomial with the algorithms from Theorem 
6.6 below. Namely, since we know, that for an integer root, say a, of the Bernstein-Sato 
polynomial of a polynomial in n > 2 variables — n + 1 < a < — 1 holds (by [33, 41]) and 
— 1 is always a root, we can run the checkRoot procedure (which is just one Grobner 
basis computation with an arbitrary ordering, see Section 6.4) starting from a = —n + 1 
to a — —2. We stop at the first affirmative answer from checkRoot or output —1 if no 
positive answer appears. Thus, one executes checkRoot at most n — 2 times. 

Algorithm 4.5 (Heuristic for AmiD n (/ a )). 

Input: / G C[xi, . . . , x n ], a G C 
Output: Ann Dn (/ Q ) 

if a G C \ (Z n [-n + 1, -1]) then 

{(<9i,...,<9„) ifa = 0, 

ker(L> n 1 ^T D n /(d 1 ,...,d n )) if a = m G N, (cf. 4.2), 
Ann Dn[s] (/ s ) | S=Q if a G (C \ Z) U (Z n (-oo, -n]), 

else (that is a G Z n [-n + 1, -1]) 
H := min{/3 G Z <0 | b f (/3) = 0} 

I Procedure 4.1 with 4.4 if a + 1 < a < — 1, 

Ami Dn (j a ) = < 

* I Procedure 4.4 and Aian Dn ^(f s ) \ s=a if — n + 1 < a < fi. 

end if 

return Ann£, n (f a ) 

4.3 Annihilator of a rational function 

In order to compute the annihilator / of a rational function - (it is known that D n /I is 
holonomic) we use the following lemma. 

Lemma 4.6. Let f,g G K[x±, . . . ,x n ] \ {0}. Consider the homomorphism of left D n - 
modules r : D n — > D r J Axm Dn {g~ l ),q i-> qf. Then Ann Dn (-) = kerr. 

Proof. For q G kerr = {q G D n \ qf G Annc^g -1 )}, (qf) • g~ x = q • (fg^ 1 ), hence 
Ann Dn (^) = kerr. D 

We compute Ann^^g -1 ) with Algorithm 4.5 above. Although in the case, when —1 
is not the minimal integer root of the Bernstein-Sato polynomial of g, we have to use 
expensive algorithms like 4.1, we know no other methods to compute the annihilator in 
Weyl algebras. Also, no general algorithm for computing a complete system of operator 

9 



equations (with operators including along partial differentiation also partial (g-)differences 
et cetera) with polynomial coefficients, annihilating a rational function, is known to us. In 
our opinion, the existence of an algorithm for Ami£) n (g _1 ) shows the intrinsic naturality 
of D-modules compared with other linear operators acting on K[x]. The algorithm is 
implemented in dmodapp.lib and the corresponding procedure is called annRat. 

Example 4.7. In this example we demonstrate the computation of annihilators of a ra- 
tional function. The procedure annRat takes as arguments polynomials in a commutative 
ring and returns a Weyl algebra (of type ring) together with an object of type ideal 
called LD (cf. Example 3.2). Note, that LD is given in a Grobner basis. 

LIB "dmodapp.lib"; 

ring r = 0, (x,y) ,dp; 

poly g = 2*x*y; poly f = x~2 - y~3; // we will compute Ann(g/f) 

option(redSB) ; option(redTail) ; // get reduced minimal GB 

def B = annRat(g,f); setring B; 

LD; // Groebner basis of Ann(g/f) 

==> LD[l]=3*x*Dx+2*y*Dy+l 

==> LD [2] =y~3*Dy~2-x~2*Dy~2+6*y~2*Dy+6*y 

==> LD [3] =9*y~2*Dx~2*Dy-4*y*Dy~3+27*y*Dx~2+2*Dy~2 

==> LD [4] =y~4*Dy-x~2*y*Dy+2*y~3+x~2 

==> LD [5] =9*y~3*Dx~2-4*y~2*Dy~2+10*y*Dy-10 

5 6-function with respect to weights for an ideal 

Let O^wG M> and consider the \/- filtration V = {V m \ m G Z} on D n with respect to w, 
where V m is spanned by {x a d@ \ —wa + w(3 < m} over K. That is, %i and di get weights 
— Wi and Wi respectively. Note that then the relation diXi = Xidi + 1 is homogeneous of 
degree 0. It is known that the associated graded ring © meZ V m /V m -\ is isomorphic to 
D n , which allows us to identify it with the Weyl algebra. 

From now on we assume, that / C D n is an ideal such that D n / 1 is a holonomic 
module. Since holonomic .D-modules are cyclic (e. g. [10]), for each holonomic D-module 
M there exists an ideal Im such that M = D n / Im as D-modules. 

Definition 5.1. Let ^ w G M> . For a non-zero polynomial 

p = 2_, c apx a dP G D n with all but finitely many c a p = 

we put m = max a ^{- wa + w(5 \ c a p ^ 0} G R and define the initial form of p with 
respect to the weight w as follows: 

in ( _ U)iTO )(p) := ^ c a/3 x a dP. 

o,/3€Nq: —wa+w/3=m 

For the zero polynomial, we set in(_ WiW )(0) := 0. Additionally, the ideal in(_ u , jW )(7) := 
K ■ {in(_ wu ,)(p) | p G /} is called the initial ideal of I with respect to w. 
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Definition 5.2. Let ^ w G IR> and s := Y^i=i w i x i^i- Then in(_ ww )(/) fl K[s] is a 
principal ideal in K[s]. Its monic generator bi tW (s) is called the global b-function of i" with 
respect to the weight w. 

Theorem 5.3. The global b-function is nonzero. 

We will give a proof of this well-known result in Section 5.2. 

Following its definition, the computation of the global 6-function of / with respect 
to w can be done in two steps: 

1. Compute the initial ideal I' of I with respect to w. 

2. Compute the intersection of I' with the subalgebra K[s\. 

We will discuss both steps separately, starting with the initial ideal. It is important 
to mention, that although this procedure has been described in [34], this approach was 
completely treated by Noro in [30], accompanied with a very impressive implementation 
in Risa/Asir. 

5.1 Computing the initial ideal 

In order to compute the initial ideal, the method of weighted homogenization is proposed 
in [30], which we will describe below. 

Let m,d6 R> . The G-algebra D\K := K(xi, . . . , x n , d\, . . . , d n , h | {xjXi = XiXj, djdi 
= didj,Xih = hxi, dih = hdi, djXi = Xidj + Si ! jh Ui+Vj }) is called the n-th weighted homoge- 
nized Weyl algebra with weights u, v, i. e. Xi and <9j get weights Ui and Vi respectively. 

For p = Y^ a s c ol j3X OL d 13 G D n we define the weighted homogenization of p as follows: 

H M (p) = J2 c ^h dcs ^ {p) - {ua+v ^x a d fS . 

This definition naturally extends to a set of polynomials. Here, deg, UiV \(p) denotes the 
weighted total degree of p with respect to weights u, v for x, d and weight 1 for h. 

For a monomial ordering -< on D n , which is not necessarily a well-ordering, we define 
an associated homogenized global ordering -<^ in D (uv) by setting h -<^ x,i,h -<^ 9, 
for all i and, 

p^ h) q if deg iuv) (p) < deg (uv) (q) 

or deg(« )W )(p) = deg (tt>t)) (g) and p lh=1 -< q lh=1 . 

Note that for u — v — (l,...,l) this is exactly the standard homogenization as in [34] 
and [9] . Analogue statements of the following two theorems can be found in [34] and [30] 
respectively. 

Theorem 5.4. Let F be a finite subset of D n and -< a global ordering. If G^ is a 
Grobner basis of (Hi u , v ){F)) with respect to -<( h \ then G^ h \ h=1 is a Grobner basis of (F) 
with respect to -<. 
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Theorem 5.5. Let -< be a global monomial ordering on D n and -<f- WjW ) the non-global 
ordering defined by 

x a d 13 -<(-w :W ) x~ f d 5 if —wa + w/3 < —wy + w5 

or —wa + w/3 = —wy + w5 and x a d 13 -< x^'d . 

lj Q( h ) i s a Grobner basis of H^ U ^{I) with respect to -<[_ \, then the set {vn^ w ^ w ^{g) \ 
g G G^} is a Grobner basis o/m(_ WiW) o)(-^(M^)(^)) with respect to -<( h \ 

Proof. Let /' G in-(-w >w ,o)(H(u,v)(I)) be (— w,w, 0)-homogeneous. There exist / G if( Uj „)(J) 
and g G G^ such that /' = in^ WiWj0 )(f) and lm^(h) (g) | lm^/o (/). Since /, g are 

(u, v )- homogeneous, we have 



y — w,w) ( — w,w) 



lm (^ (g) = ha^( h )(m { - WiWfi )(g)), lm w (f) = lm^ w (in ( _ W)Wi0 )(/)), 

which concludes the proof. D 

Summarizing the results from this section, we obtain the following algorithm to com- 
pute the initial ideal. 

Algorithm 5.6 (Initialldeal). Input: / C D n a holonomic ideal, O^wG ^>o> ~^ a 

global ordering on D n , u,v G IR> 
Output: A Grobner basis G of in(_ wu ,)(J) with respect to -< 

-O v= the homogenized ordering as defined in theorem 5.5 

( — W,W) ° o 

G^ h > := a Grobner basis of H( U>V )(I) with respect to <\_ ww \ 
return G = m^ w>w )(G w \ h=1 ) 

5.2 Intersecting an ideal with a principal subalgebra 

We will now consider a much more general setting than needed to compute the global 
6-function. Let A be an associative IK-algebra. We are interested in computing the 
intersection of a left ideal J G A with the subalgebra K[s] of A where s G A is an arbitrary 
non-constant element. This intersection is always generated by one element since K[s] is 
a principal ideal domain. In other words, we want to find the monic generator b G A such 
that (b) = JflK[s]. 

For this section, we will assume that there is an ordering on A such that there exists 
a finite left Grobner basis G of J. 

Then we can distinguish between the following four situations: 

1. No leading monomials of elements in G divide the leading monomial of any power 
of s. 

2. There is an element in G whose leading monomial divides the leading monomial of 
some power of s. In this situation, we have the following sub-situations. 
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2.1. J ■ s C J and dim K (EndA(A/</)) < oo. 

2.2. One of the two conditions in 2.1. does not hold. 

2.2.1. The intersection is zero. 

2.2.2. The intersection is not zero. 

Lemma 5.7. If there exists no g G G such that lm(g) divides lm(s fe ) for some k G N , 
then Jf]K[s] = {0}. 

The lemma covers the first case above. In the second case however, we cannot in 
general state whether the intersection is trivial or not as the following example illustrates. 

Remark 5.8. The converse of the previous lemma is wrong. For instance, consider K[x, y] 
and J = (y 2 + x). Then J fl K.[y] = {0} while {y 2 + x} is a Grobner basis of J for any 
ordering. 

In situation 2.1. though, the intersection is not zero as the following lemma shows, 
inspired by the sketch of the proof of Theorem 5.3 in [34]. 

Lemma 5.9. Let J ■ s C J and dim K (EncU(A/J)) < °°. Then J H K[s] ^ {0}. 

Proof. Consider the right multiplication with s as a map A/ J — > A/ J which is a well- 
defined A-module endomorphism of A/ J as a — a' G J implies that (a — a')s G J ■ s C J, 
which holds by assumption for all a, a' G A. Since End^A/J) is finite dimensional, 
linear algebra guarantees that this endomorphism has a well-defined non-zero minimal 
polynomial \x. Moreover, \x is precisely the monic generator of J fl K[s] as /x(s) = [0] in 
A/ J, hence /x(s) G J fl K[s], and deg(yu) is minimal by definition. □ 

Remark 5.10. In particular, the lemma holds if A/ J itself is a finite dimensional A- 
module. In the case where A is a Weyl algebra and A/ J is a holonomic module, we know 
that dim K (End A (A/J)) is finite (cf. [34]). 

For situation 2.1., we have reduced our problem of intersecting an ideal with a sub- 
algebra generated by one element to a problem from linear algebra by the proof of the 
lemma, namely to the one of finding the minimal polynomial of an endomorphism. 

Proof of Theorem 5.3. Let O^wG ^>o> J := m (-«>,w)(-0 f° r a holonomic ideal / C D n 
and s := Y^i=i w i x i^i- Without loss of generality let ^ p = Y^ a B c a,P xa ® G J be 
(—w, u>)- homogeneous. Then we obtain for every monomial in p by using the Leibniz rule 

x a d p x i d i = x a+ei d? +ei + (3 t x a dP = (d^ +1 - {a t + l)xf)^-d p + ^x a d^ 

= (dpi - (a, + 1) + A)x^ = (x t d t -a t + f3 t )x a d?. 

Put m = —wa + w/3 for some term c a ^x a d^ in p where c Q;/ g is non-zero. Since p is 
(— 117,1/7)- homogeneous, m does not depend on the choice of this term. Hence, 

n n 

p. s = p^2 WiXidi = ^2 w i ^2( X A -(Xi + I3i)c a ^x a d p 

i=l i=l a, /3 

n 

= S -p + ^2^2 w i{~ a i + Pi)c a ,/3 xa dP = (s + m) ■ p G J. 

i=l a,{3 
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Since D n /J is holonomic (cf. [34]) and J ■ s C J, Remark 5.10 and Lemma 5.9 yield the 
claim. □ 

If one knows in advance that the intersection is not zero, the following algorithm 
terminates. 

Algorithm 5.11 (principallntersect). 

Input: s E A, J C A a left ideal such that J n K[s] ^ {0}. 
Output: b E K[s] monic such that J D K[s] = (b) 

G := a finite left Grobner basis of J (assume it exists) 

i := 1 

loop 



if there exist a , . . . , Oj-i G K such that NF(s\ G) + Yl]=o a j NF(s J ', G) = then 



return 6 := s* + X^=o a i s ' ? 
else 

i :=i + l 
end if 
end loop 

Note that because NF(s J , G) + Yll^o a i NF(s J , G) = is equivalent to s l + Y^jlo a j s ^ *= 
J, the algorithm searches for a monic polynomial in K[s] that also lies in J. This is done 
by going degree by degree through the powers of s until there is a linear dependency. 
This approach also ensures the minimality of the degree of the output. The algorithm 
terminates if and only if J D K[s] ^ {0}. Note that this approach works over any field. 

The check whether there is a linear dependency over IK between the computed normal 
forms of the powers of s is done by the procedure linReduce in our implementation. 

5.2.1 An enhanced computation of normal forms 

When computing normal forms of the form NF(V, J) like in algorithm 5.11 we can speed 
up the reduction process by making use of the previously computed normal forms. 



Lemma 5.12. Let A be a K.-algebra, J C A a left ideal and let f G A. For iGN put 

lc(ri<ji) 



Ti = NF(/ J , J), qi = f l — Ti G J and a = , ( q -, provided r\qi ^ 0. For r\qi = we put 



Ci = 0. Then we have for alii EN 

r m = NF(/r<, J) = NF([f - r u n] Ci + nn, J). 

As a consequence, we obtain the following result for some K-algebras of special im- 
portance. 

Corollary 5.13. If A is a G-algebra of Lie type (e. g. a Weyl algebra), then 

r i+1 = NF(/ n , J) = NF([f - n, n] + Ti r u J) holds. 
If A is commutative, we have r\ + i = NF(rjri, J) = NF(n, J) %+1 = NF(r^ +1 , J). 

Note, that computing Lie bracket [/, g] both in theory and in practice is easier and 
faster, than to compute [f,g]&sf-g — g-f, see e. g. [22]. 
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5.2.2 Applications 

Apart from computing global 5-functions, there are various other applications of Algo- 
rithm 5.11. 

Solving Zero-dimensional Systems. Recall that an ideal / C K.[x±, . . . ,x n ] is called 
zero- dimensional if one of the following equivalent conditions holds: 

• K.[xi, . . . , x n ]/I is finite dimensional as a K- vector space. 

• For each 1 < i < n there exist ^ fi G / fl K[xj]. 

• The cardinality of the zero-set of / is finite. 

In order to compute the zero-set of J, one can use the classical triangularization 
algorithms. These algorithms require to compute a Grobner basis with respect to some 
elimination ordering (like lexicographic one), which might be very hard. 

By Algorithm 5.11, a generator of / fl K[xi\ can be computed without these expensive 
orderings. Instead, any ordering, hence a better suited one, may be freely chosen. 

A similar approach is used in the celebrated FGLM algorithm (cf. [12]). 

Computing Central Characters and Algebraic Dependence. Let A be an asso- 
ciative K-algebra. Intersection of a left ideal with the center of A, which is isomorphic to 
a commutative ring, is important for many algorithms, among other for the computation 
of central character decomposition of a finitely presented module (cf. [19] for the theory 
and [1] for an example with Principal Intersection). In the situation, where the center 
of A is generated by one element (which is not seldom), we can apply Algorithm 5.11 
to compute the intersection (known to be often quite nontrivial) without engaging much 
more expensive Grobner basis computation, which use elimination. 

Example 5.14. Consider the quantum algebra U'Jsos) (as defined by Fairlie and Odesskii) 
for q 2 being the n-th root of unity. It is known, that then, in addition to the single gen- 
erator C of the center present over any field, three new elements Zj, depending on n will 
appear. Since ^(503) has Gel'fand-Kirillov dimension 3, four commuting elements in it 
obey a single polynomial algebraic dependency (the ideal of dependencies in principal). 
Computing such a dependency is a very tough challenge for Grobner bases. But as we 
see, it is quite natural to apply Principal Intersection. 

LIB "ncalg.lib"; LIB "bfun.lib"; 

def A = makeQso3(5) ; // below Q~2 is the 5th root of unity 

setring A; 

// central elements, depending on Q in their classical form: 

ideal I = x5+(Q3-Q2+2) *x3+(Q3-Q2+l) *x, y5+(Q3-Q2+2) *y3+(Q3-Q2+l) *y , 

z5+ (Q3-Q2+2) *z3+ (Q3-Q2+1) *z ; 
I = twostd(I) ; // two-sided Groebner basis 

poly C = 5*xyz+(4q3-3q2+2q-l)*x2+(-q3+2q2-3q+4)*y2+(4Q3-3Q2+2Q-l)*z2; 
poly v = vec2poly (plntersect (C, I) , 1) ; // present vector as poly 
poly t = subst (v,x,C) ; t; // t as a polynomial in C of size 42 
==> 3125*x5y5z5+(3125q3-3125q2+6250)*x5y5z3+(3125q3-3125q2+6250)*x5y3z5+ . . . 
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// present matrix of cof actors of t as an element of I : 

matrix T = lift(I,t); 

poly a = 125*(25*I[l]*I[2]*I[3] + (Q3-7Q2+8Q-4)*(I[l]~2+I[2]-2+I[3]-2)); 

a-t ; // a expresses t in the subalgebra gen. by I[1..3] 

==> 

// define univariate ring over algebraic extension: 

ring r = (0,Q),c,dp; minpoly = Q4-Q3+Q2-Q+1 ; 

poly v = fetch(A,v); // map v from A to a univariate poly in c 

f actorize(v) ; 

The latter factorization delivers the final touch to the answer: the algebraic depen- 
dency is described by the equation C 2 ■ (C + 4g 3 — 3q 2 — 3q + 4) • (C + 3g 3 — q 2 — q + 3) 2 = 
3125 • Z t Z 2 Z 3 + 125 • (q 3 - 7q 2 + 8q - 4) • [Z 2 + Z 2 + Z 2 ). 

6 Bernstein-Sato polynomial of / 

6.1 Global Bernstein-Sato polynomial 

One possibility to define the Bernstein-Sato polynomial of a polynomial / G K[xi, . . . , x n ] 
is to apply the global 6-function for specific weights. 

Definition 6.1. Let bj w (s) denote the global 6-function of the univariate Malgrange 
ideal If of / (cf. Section 3) with respect to the weight vector w = (1, 0, . . . , 0) G IR n+1 , 
that is the weight of d t is 1. Then bf(s) = (—1) If ' w bj f:W (—s — 1) is called the global 
b-function (Bernstein- Sato polynomial) of /. 

By Theorem 5.3, bf(s) ^ holds. Moreover, it is known that all roots of 6/(s) 
are negative rational numbers. Kashiwara proved this result for local Bernstein-Sato 
polynomials over C [18]. This fact together with Theorem 6.3 below and classical flatness 
properties imply the claim for the global case over an arbitrary field of characteristic 0. 

The following theorem gives us another option to define the Bernstein-Sato polynomial. 

Theorem 6.2 ([4], see also [34, Lemma 5.3.11]). The Bernstein- Sato polynomial bf(s) of 
f is the unique monic polynomial of minimal degree in K[s] satisfying the identity 

P • / s+1 = bf(s) ■ f s for some operator P G -D„[s]. 

Since P ■ f — bf(s) G AnnD n [ s ](f s ) holds, 6/(s) is the monic polynomial satisfying 

(b f ( S )) = Ann Dn[s] (f) + (f)nK[s}. (1) 

Summarizing, there are several choices for computing the Bernstein-Sato polynomial: 

1. Compute either 

(a) J = in(_ W)tu )(i/) or 
(b)J = Ann Dn[a] (f') + {f). 

2. Intersect J with K[s] by 

(a) the classical elimination-driven approach or 

(b) using Algorithm 5.11. 
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It is very interesting to investigate the approach for the computation of Bernstein-Sato 
polynomial, which arises as the combination of the two methods: 

1. Anii£) n [ s ](/ S ) via Briangon-Maisonobe (cf. [20]), 

2. (Ann( Dn[s] f s ) + (/)) n K[s] via Algorithm 5.11. 

For an efficient computation of mt_ w>w -j(If) using the method of weighted homog- 
enization as described in Section 5.1, Noro proposes [30] to choose the weights u = 
(deg u (/),«i,...,u„), v = (l,deg u (/)-tii + l,...,deg tt (/)-u n + l), such that the weight 
of t is deg u (/) and the weight of d t is 1. Here, u G M™ is an arbitrary vector and deg u (/) 
denotes the weighted total degree of / with respect to u. The vector u may be chosen 
heuristically in accordance to the shape of / or by default, one can set u — (1, . . . , 1). 

6.2 Implementation 

For the computation of Bernstein-Sato polynomials, we offer the following procedures in 
the Singular library bfun.lib: 

bfct computes mr_ WjW <\(If) using weighted homogenization with weights u, v for an 
optional weight vector u (by default u = (1, . . . , 1)) as described above, and then uses Al- 
gorithm 5.11, where the occurring systems of linear equations are solved by the procedure 
linReduce. 

bfct Ann computes Annr> n [ s ](f s ) via Briancon-Maisonobe and intersects Ann£> n [ s ]{f s ) + 
(/) with K[s] analogously to bfct. 

bf ctOneGB computes the initial ideal and the intersection at once using a homogenized 
elimination ordering, a similar approach has been used in [16]. 

For the global 6-function of an ideal / C D n , bf ctldeal computes m^ w>w )(I) using 
standard homogenization, i. e. weighted homogenization where all weights are equal to 1, 
and then proceeds the same way as bfct. Recall that D n /I must be holonomic as in [34]. 

All these procedures work as the following example illustrates for bfct and the hyper- 
plane arrangement xyz(z — y)(y + z). 

LIB "bfun.lib"; 

ring r = 0, (x,y,z) ,dp; // commutative ring 

poly f = x*y*z*(z-y)*(y+z) ; 

list L = bfct(f) ; 

print (matrix (L [1] )) ; // the roots of the BS-polynomial 

==> -1,-5/4,-3/4,-3/2,-1/2 

L [2] ; // the multiplicities of the roots above 

==> 3,1,1,1,1 

6.3 Local Bernstein-Sato Polynomial 

Here we are interested in what kind of information one can obtain from the local b- 
functions for computing the global one and conversely. In order to avoid theoretical 
problems we will assume in this paragraph that the ground field K = C. 

Several algorithms to obtain the local 6-function of a hypersurface / have been known 
without any Grobner bases computation but under some conditions on /. For instance, 
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it was shown by Malgrange [24] that the minimal polynomial of —d{t acting on some 
vector space of finite dimension coincides with the reduced (local) Bernstein polynomial, 
assuming that the singularity is isolated. 

The algorithms of Oaku [32] used Grobner bases for the first time. Recently, Nakayama 
presented some algorithms, which use the global 6-function as a bound and obtain a local 
6-function by Mora resp. approximate division [27], see also the work of Nishiyama and 
Noro [29]. 

Theorem 6.3 (Briangon-Maisonobe (unpublished), Mebkhout-Narvaez [26]). Letbf j p(s) 
be the local b-function of f at the point PeC" and bf(s) the global one. Then it is ver- 
ified that bf(s) = lcmpgo b ftP (s) = lcm PeS(/ ) b fjP (s), where £(/) = V({f, J£, . . . , J£)) 
denotes the singular locus ofV(f). 

Remark 6.4. Assume, that S(/) consists of finitely many isolated singular points (the 
dimension of the corresponding defining ideal is 0). Then the computation of the global 
6-function with Theorem 6.3 becomes effective. Moreover, one needs just an algorithm 
for computing the local 6-function of a hypersurface, having an isolated singularity at the 
origin. 

The Singular library gmssing.lib, developed and implemented by M. Schulze [36], 
contains the procedure bernstein, which computes the local 6-function at the origin. It 
returns the list of roots and corresponding multiplicities. 

Example 6.5. Let C be the curve in C 2 given by / = (x 3 — y 2 )(3x — 2y — l)(x + 2y). 
This curve has three isolated singular points pi = (0, 0), p 2 = (1, 1) and p 3 = (1/4, —1/8). 

LIB "gmssing.lib"; 

// note that one must use a local ordering for calling 'bernstein' 

ring r = 0,(x,y),ds; // ds stands for a local degrevlex ordering 

poly f = (x~3-y~2)*(3x-2y-l)*(x+2y) ; 

list L = bernstein(f ) ; // local b-function at the origin p_l 

print (matrix (L [1] ) ) ; 

==> -11/8,-9/8,-1,-7/8,-5/8 

L[2]j 

==> 1,1,2,1,1 

Moving to the corresponding points we also compute &/,p 2 (s) and 6/ : p 3 (s). 

bf,Pi(s) = ( s + l ) 2 ( s + 5 / 8 )( s + 7 /8)(s + 9/8) (s + 11/8) 
W*) = (s + l) 2 (s + 3/4)( S + 5/4) 
W*) = (s + l) 2 (s + 2/3)( S + 4/3) 

From this information and using Theorem 6.3, the global 6-function is 

bf(s) = (s + 2/3)(s + 5/8)(s + 3/4)(s + 7/8)(s + l) 2 (s + 4/3)(s + 5/4)(s + 9/8)(s + 11/8). 

Moreover, gmssing. lib, allows one to compute invariants related to the Gauss-Manin 
system of an isolated hypersurface singularity. 
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In the non-isolated case the situation is more complicated. For computing the local b- 
function in this case (which is important on its own) we suggest using two methods: Take 
the global ^-function as an upper bound and a local version of the checkRoot algorithm, 
see below. Another method is to use a local version of principallntersect, which is 
under development. Despite the existence of many algorithms, the effectiveness of the 
computation of local 6-functions is still to be drastically enhanced. 

6.4 Partial knowledge of Bernstein-Sato polynomial 

As we have mentioned, several algorithms for computing the 6-function associated with 
a polynomial have been known. However, in general it is very hard from computational 
point of view to obtain this polynomial, and in the actual computation a limited number 
of examples can be treated. For some applications only the integral roots of b/(s) are 
needed and that is why we are interested in obtaining just a part of the Bernstein-Sato 
polynomial. 

Recall the algorithm checkRoot for checking whether a rational number is a root of 
the 6-function of a hypersurface from [20]. Equation (1) was used to prove the following 
result. 

Theorem 6.6. ([20]) Let R be a ring whose center contains K[s] as a subring. Let 
us consider q(s) G K[s] a polynomial in one variable and I a left ideal in R satisfying 
I n K[s] ^ 0. Then (I + R(q(s))) D K[s] = I D K[s] + K[s](g(s)). In particular, using the 
above equation (1), we have 

( Aim Dn[s] (f s ) + D n [s] ■ (f,q(s))) nK[s] = (b f (s),q(s)). 

As a consequence, let m a be the multiplicity of a as a root ofbf(-s) and let us consider 
the ideals Ji = Anii£) n [ s ](/ S ) + (/, (s + a) t+1 ) C -D n [s], i = 0, . . . ,n, then [m a > i •<=>- 
{s + aYiJi}. 

Once we know a system of generators of the annihilator of f s in -D n [s], the last theorem 
provides an algorithm for checking whether a given rational number is a root of the b- 
function of / and for computing its multiplicity, using Grobner bases for differential 
operators. 

This algorithm is much faster, than the computation of the whole Bernstein polynomial 
via Grobner bases, because no elimination ordering is needed for computing a Grobner 
basis of Ji, once one knows a system of generators of Ama.D n [ s ](f s ). Also, the element 
(s + a) t+1 , added as a generator, seems to simplify tremendously such a computation. 
Actually, when % — it is possible to eliminate the variable s in advance and we can 
perform the whole computation in D n . Let us see an example. 

Example 6.7. Let A be the matrix given by 

X\ %2 X3 X4 

A = ( £5 xq xj x 8 

Xg £10 £11 X\2 
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Let us denote by A«, % — 1, 2, 3, 4, the determinant of the minor resulting from deleting the 
i-th column of A, and consider / = Ai A2A3A4. The polynomial / defines a non-isolated 
hypersurface in C 12 . Therefore, from [33] (see also [41]), the set of all possible integral roots 
of bf(-s) is {11,10,9,8,7,6,5,4,3,2,1}. It is known that Ann DnH (/ s ) = Ann^ [s] (/ S ) 
(see Section 8) and this fact can be used to simplify the computation of the annihilator. 

LIB "dmod.lib"; 

ring R = 0, (xl ,x2,x3,x4,x5,x6,x7,x8,x9,xl0,xll ,xl2) ,dp; 

matrix A [3] [4] = xl,x2,x3,x4, x5,x6,x7,x8, x9,xl0,xll ,xl2; 

poly Deltal = det (submat (A, 1 . .3,intvec(2,3,4) ) ) ; 

... // analogous for Delta2 . . . Delta4 

poly f = Deltal*Delta2*Delta3*Delta4; 

def D = Sannf slog(f ) ; setring D; // logarithmic annihilator 

poly f = imap(R,f); number alpha = 11; 

checkRootl(LDl,f .alpha) ; 

==> 

Using the algorithm checkRoot we have proved that the minimal integral root of bf(s) 
is — 1. This example was suggested by F. Castro- Jimenez and J. M. Ucha for testing the 
Logarithmic Comparison Theorem. A nice introduction to this topic can be found, for 
instance, in [39]. 

Let g be the polynomial resulting from / by substituting xi,x 2 ,x 3 ,x 4: ,x 5 ,Xq with 1. 
One can show that b g (s) divides bf(s) (see [21] for details). Using the checkRoot algorithm 
we have found that (s + l) 4 (s + 1/2) (s + 3/2) (s + 3/4) (s + 5/4) is a factor of b g (s) and 
therefore a factor of bf(s). 

Remark 6.8. Using the notation from Section 5, given a holonomic D-module D/I, 
it is verified that (in^ W;W )(I) + (q(s))) fl K[s] = (bi jW (s),q(s)), although Theorem 6.6 
cannot be applied, since s = Y2 w i x idi does not commute with all operators. For some 
applications like integration and restriction the maximal and the minimal integral root 
of the b- function of / with respect to some weight vector have to be computed, see [34]. 
However, the above formula cannot be used to find the set of all integral roots, since no 
upper/lower bound exists in advance. For instance, as it was suggested by N. Takayama, 
/ = (tdt + k) , k G Z is .Di-holonomic and one has in(_ 11 )(J) DC [5] = (s + k) with s = tdt- 

We close this section by mentioning that there exist some well-known methods to 
obtain an upper bound for the Bernstein-Sato polynomial of a hypersurface singularity 
once we know, for instance, an embedded resolution of such singularity [18]. Therefore 
using this result by Kashiwara and the checkRoot algorithm, it is possible to compute 
the whole Bernstein-Sato polynomial without elimination orderings, see Example 1 in 
[20]. We investigate different methods in conjunction with the further development of the 
checkRoot family of algorithms in [21]. 

7 Bernstein operator of / 

We define the Bernstein-Sato polynomial bf(s) to be the monic generator of a principal 
ideal, hence it is unique. But the so-called 5-operator P(s) G D n [s] from Theorem 6.2 is 
not unique. 
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Proposition 7.1. Let G be a left Grobner basis o/Annj) n [ s ](/ s+1 ) and define a Bernstein 
operator to be the result of the reduced normal form NF(P(s),G) of some B-operator 
P(s). Then, for a fixed monomial ordering on D n [s], the Bernstein operator is uniquely 
determined. 

Proof. Suppose that there is another Q(s) G -D n [s], such that the identities 

P(s)f s+1 = b f (s)f s and Q(s)f s+1 = b f (s)f s hold in the module D n [s}/ Ann Dn[s] (f). 

Then (P(s) - Q(s))f s+1 = 0, that is P(s) - Q(s) e Ann Dn[s] (f s+1 ). Hence, the set 
{R(s) G D n [s] | R(s)f s+1 = bf(s)f s } can be viewed as an equivalence class and we 
can take the reduced normal form of any such operator (with respect to G) to be the 
canonical representative of the class. Since the reduced normal form with respect to a 
fixed monomial ordering on D n [s] is unique, so is the Bernstein operator NF(i?(s), G). □ 

Note, that we can obtain the left Grobner basis of Ann^ r s i(/ s+1 ) via substituting s 
with s + 1 in the left Grobner basis of Anno n [ s ](f s )- 

One can compute the Bernstein operator from the knowledge of AmiD n [ s ](/ s ) and 6/(s) 
by the following methods. 

7.1 5-operator via lifting 

The algorithm Lift(F, G) computes the transformation matrix, expressing the set of 
polynomials G via the set F, provided (G) C (F). It is a classical application of Grobner 
bases. 

Lemma 7.2. Suppose, that Awn.D n [ s ](f s ) is generated by hi, ... , h m and bf(s) is known. 
The output of Lift({/, hi, . . . , h m }, {&/($)}) is the 1 x (m + 1) matrix (a, bi, . . . , b m ). 
Then a B-operator is computed as P(s) = NF(a, Anno n [ s ]{f s )). 

Proof. Because of Equation (1), if (a, bi, . . . , b m ) is the output of Lift as in the statement, 

m 

bf(s) = af + \^ bihi holds, 

i=l 

hence the first element of such a matrix is a -B-operator. Thus, the Bernstein operator is 
obtained via NF(a, Ann Dn [ s ](/ s+1 )). □ 

However, we have to mention, that the Lift procedure is quite expensive in general. 
Note that another method for the computation of a 5-operator using lifting techniques 
is given by applying Algorithm 8 of [29] with a(x) = 1. 

7.2 5-operator via kernel of module homomorphism 

1. Consider the -Cv[s]- module homomorphism 

y? : D n [s] — ► D n [ S }/(Aim Dn[s] (f) + (b f (s))), 1 h- /, 
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then for u G kenp, uf G Ann Dn ^(f s ) + (bf(s)). That is, there exist a, bi G -D n [s], such 
that 



m 



uf = abf(s) + 2_^bihi, 

i=l 

However, we are interested in such u, that a G K. This is possible, but the 2nd method 
above proposes a more elegant solution. Also one has to say, that in this case we have 
to compute a Grobner basis of AnnD„[ s ](/ s ) + (bf(s)) as an intermediate step and also 
the kernel of a module homomorphism with respect to the latter. This combination is, 
in general, quite nontrivial to compute. In the Grobner basis computation a monomial 
ordering, preferring x, dx over s seems to be better because of numerous applications of 
the Product Criterion. 

2. Consider the Z) n [s]-module homomorphism 

ti : D n [s} 2 — ► D n [s\/ Ann DnH (/ s ), ei ^ b f (s),e 2 m- /. 

Then ker$ = {(u,v) T G -D n [s] 2 | ubf(s) +vfE Aim.D n [ s ](f s )} is a submodule of D n [s} 2 . 
Indeed, ker $ has many generators. In order to get a vector of the form (k, u(s)) for k G K, 
we perform another Grobner basis computation for a submodule with respect to a module 
monomial ordering, giving preference to the first component over the second one. Since 
in the reduced basis there is a single element of the form (k,v(s)) C ker$ with k ^ 0, it 
follows that P(s) = v (s)k' 1 . 

This algorithm is implemented in dmod.lib as operatorModulo. The approach via 
lifting is used in the procedure operatorBM, which computes all the Bernstein data. The 
procedures can be used as follows. 

Example 7.3. Consider the Reiffen curve / = x 2 + y 3 + xy 2 G K[x, y]. At first we use 
operatorBM and compare the length of an 5-operator computed via lift to the length of 
the Bernstein operator. 

LIB "dmod.lib"; 

ring r = 0, (x,y) ,dp; 

poly F = x~2 + y~3 + x*y~2; 

def A = operatorBM (F) ; setring A; 

size (PS); // size of B-operator 

==> 238 

ideal LD2 = subst (LD,s ,s+l) ; // Ann(F~{s+l}) 

LD2 = groebner (LD2) ; // LD is not a Groebner basis 

poly PS2 = NF(PS,LD2); size(PS2); // size of Bernstein operator 

==> 41 

So, computing with lifting potentially computes much longer operators. Let us compare 
with operatorModulo. 

poly PS3 = operatorModulo(F,LD,bs) ; size(PS3); 
==> 41 

The size of the operator, returned by operatorModulo need not be minimal (e. g. by dis- 
abling some of the interactive options of Singular one can get a polynomial of length 50 
in this example), but it is in general much shorter, than the one, delivered by operatorBM. 
Let us check the main property of the .B-operator and print its highest terms: 
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NF(PS3*F - bs, subst(LD2,s,s-l)); 

==> 

108*PS2; // i.e. the Bernstein operator 

==> 6*x*Dx~2*Dy+9*y*Dx~2*Dy-2*x*Dx*Dy~2-4*y*Dx*Dy~2-y*Dy~3+ . . . 

In the last line we see the terms of highest degree with respect to dx, dy. 

7.3 Grobner free method 

As a consequence of Theorem 6.2, one obtains that P ■ f — bf(s) G Ann£> n [ s ]{f s ) and 
P,bf(s) <£ Ann Dn [ s ](/ s ). If we fix an ordering such that &/(s) = NF(6/(s), Ann Dn ^(f s )) 
holds, we may rewrite this relation to &/(s) = NF(P • /, AnnD n [ s ](/ s )). Hence, we can 
compute P by searching for a linear combination of monomials m G D n [s] that satisfy this 
equality when multiplied with / from the right side. Using the results from the beginning 
of this section, one only needs to consider monomials which span D n / Aim.D n [ s ]{f s+1 ) as 
K-vector space. We get the following algorithm. 

Algorithm 7.4. 

Input: / G K[xi, . . . ,x n ], the Bernstein-Sato polynomial 6/(s) of / 
Output: P G -D n [s], the Bernstein operator of / 
d:=0 
loop 

Md := {m G D n [s] \ m monomial, deg(m) < d, fm(p) \m\/p G Ann^ r s i(/ s+1 )} 
if there exist a m G K such that bf(s) = J2 meM a m NF(m • /, Ann Dn ^(f s )) then 
return P := Yl a m m — bf(s) 

m£M d 

else 

d:= d+ 1 
end if 
end loop 

The search for the coefficients a m can be done using linReduce (cf. Algorithm 5.11) 
as one is in fact looking for a linear dependency between the Bernstein-Sato polynomial 
and the elements m • f in the vector space D n [s]/ Ann^ m(/ s ). 

Remark 7.5. Note, that Algorithm 7.4 can be extended to one searching for both B- 
operator and Bernstein-Sato polynomial simultaneously. We have to mention, that both 
algorithms of this kind are well suited for the search of operators and Bernstein-Sato 
polynomials in the case, when both of them are of relatively low total degree. 

7.4 Computing integrals and zeta functions 

Given a simplex C C K™ (for K. — M, C) and / G K[x±, . . . ,x n ], we can define ((s) : = 
j c f(x) s dx. Since P(s) • / s+1 = bf(s)f s , we obtain 

C(s) = J f(x) s dx = ^y / P(s) • f(xy +1 dx 
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expanding the latter with e. g. the chain rule, we come to an in general inhomogeneous 
recurrence relation for C( s )> which involves coefficients in K[s]. Since P(s) is globally 
defined (and is, of course, independent on C), one can obtain a generic formula for all 
integrals of this type. 

Example 7.6. Let / = x 2 — x G K[x]. Then the Bernstein operator reads as P(s) = 
(2x — l)d x — 4(s + 1) and 6/(s) = s + 1. Any simplex in IK 1 is the interval [a, b] =: C. 

C(s) = / f(x) s dx = -j- / ((2s - l)3 x - A(s + 1)) . f(x) s+l dx 
Jc s + 1 J c 

1 f (2x - l)d x • f(x) s+1 dx - A((s 4 1 ; 



l-/c 



By the chain rule, J c (2x - l){d x • f(x) s+1 )dx - (2x - l)/(a;) s+1 \c -2 J c f(x) s+1 )dx, 
hence 



1 .... x.+i. 2 



^ s ) = 7TI • ^ ~ l )f( x ) \c -J~l^ s + !) - 4 C( S + 1). 



5 + 1 v " v ' '5 + 

and thus 

(4s + 6)C(5 + 1) + (5 + l)C(s) = (26 - 1)(6 2 - 6) s+1 - (2a - l)(a 2 - a) s+1 

The right hand side, say R(s), satisfies the homogeneous recurrence R(s + 2) — (a 2 — a + 
6 2 - 6).R(s + 1) + (a 2 - a)(6 2 - b)R(s) = of order 2. Substituting the left hand side into 
it, we obtain a homogeneous recurrence with polynomial coefficients of order 3: 

(a 2 - a)(b 2 -b)(s+ l)C(s) - ((s + 2)(a 2 - a + b 2 - b) - (As + 6)(a 2 - a)(6 2 - b))((s + 1) 
-((4s + 10)(a 2 -a + b 2 -b)-(s + 3))C(s + 2) + (4s + 14)C(s + 3) = 0. 

To guarantee the uniqueness of a solution to this equation, we need to specify 3 initial 
values, which can be easily done. However, such recurrences very seldom admit a closed 
form solution, thus most information about ((s) is contained in the recurrence itself. 

8 Logarithmic annihilator of / 

Given a polynomial / G K[x] = K.[x\, . . . , x n ], consider the left ideal Ann^\,(/ S ) C D n [s] 
generated by those differential operators P(s) G D n [s\ of total order (in the partials) less 
than or equal to one, which annihilate f s . This ideal is clearly contained in Anno n [ s ](f s ) 
and can be generated by elements of the form 

P(s) = a (x, s) + ai(x, s)d xi H h a n (x, s)d Xn G D n [s], 

where (ao, Oi, • • • , a n ) G syz K r XjS i(/, s-g^-, ■ ■ ■ , Sq£-). Therefore, for each / G K[x] one can 

compute, by using Grobner bases in K[x, s], a system of generators of Ann D , s Af s ). The 
corresponding procedure in dmod.lib is called Sannfslog. Let us see the Reiffen curve 
/ = x A + y 5 + xy 4 as an example with Singular. 
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LIB "dmod.lib"; 

ring R = 0, (x,y) ,dp; 

poly f = x~4+y~5+x*y~4; 

def A = Sannf slog(f ) ; setring A; LD1; 

==> LD1 [1] =4*x~2*Dx+5*x*Dx*y+3*x*y*Dy-16*x*s+4*y~2*Dy-20*y*s 

==> LD1 [2] =16*x*Dx*y"2-125*x*Dx*y-4*x~2*Dy+4*Dx*y"3+5*x*y*Dy+12*y~3*Dy-100*y~2*Dy 

-64*y~2*s+500*y*s 
// now we compute the whole annihilator with Sannfs and compare 
setring R; def B = Sannf s(f); setring B; 
map F = A,x,Dx,y,Dy,s; 
ideal LD1 = F(LD1) ; 
LD1 = groebner(LDl) ; 
simplify ( NF(LD,LD1), 2); 
==> _[l]=36*y~3*Dx~2-36*y~3*Dx*Dy+1125/4*x*y*Dx~2-315/4*x*y*Dx*Dy+ . . . 

And the latter polynomial is not an element of Ann^ y(/ s ) but of Ann^ u(/ s ) = 
Ann Dn[s] (/ s ). 

8.1 The annihilator up to degree k 

More generally, for a given k > 1 one can consider the left ideal Ann^A ,(/ s ) C D n [s] 
generated by the differential operators P{s) G D n [s] of total order less than or equal to 
k, such that P(s) annihilate f s . The tower of ideals 

Ann£ [s] (r) C • • • C Ann^f) = Ann Dn[s] (f) 

has been recently studied by Narvaez in [28]. It is an open problem to find the minimal 
integer fc satisfying the above condition without computing the whole annihilator. 

Computationally the annihilator up to degree k can be obtained using Grobner bases 
in K[x, s] as follows. Consider P(s) = Yl\B\<k a P^ e Annj-A, (/ s ) an< ^ ^ 9p( x -> s ) e 
K[x, s] be the polynomial defined by the formula d@ ■ f s = gp • f s ~\@\. Then (Q>p)\0\<k £ 
syz((yf / g/ fc_ l /3 l)| / g|< fc . Eventually, the polynomials gp(x,s) can be computed using the ex- 
pression given in Lemma 8.1. 

Given f3 = (ft, ...,£„) G N n , as usual \/3\ = ft + • • • + ft, ft = ft! • • • ft! and 
A^ = jnpP. A partition of /3 is a way of writing /3 as a sum of integral vectors with 
non-negative entries. Two sums which only differ in the order of their summands are 
considered to be the same partition. If /3 — 0\ + . . . + a^ with Cj 7^ 0, then a is said to 
be a partition of length k. The set of all partitions of /3 (resp. of length k) is denoted by 
V(/3) (resp. V(/3;k)). Obviously V(/3) = ujfi^ft A;). Finally, we write 1(a) := (Lr)r, 
where l UT is the number of times that r appears in a. 

Lemma 8.1. Using the above notation for all non-zero (3 G N" we have, 

W / x 1 

fc=l ^ ' a£T(0;k) ^ >' 
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This formula was suggested by Narvaez and can be proved by induction on \j3\. Similar 
expressions appear in [26, Prop. 5.3.5] and [25, Prop. 2.3.2]. 

However, despite this almost closed form, the set of polynomials, between which we 
have to compute syzygies, is growing fast and the size of polynomials increases. This 
results in quite hard computations even with the mentioned enhancements. 

9 Bernstein- Sato ideals for f = fi • ■ ■ ■ • f m 

Using the results from [13], which we confirmed through intensive testing (cf. [20]), 
it follows, that the method by Briangon-Maisonobe is the most effective one for the 
computation of s-parametric annihilators where / = f\. Because of the structure of 
annihilators in the situation / = f\ ■ ■ ■ f p , p > 1, basically the same principles stand 
behind the corresponding algorithms. Hence, we decided to implement only Briangon- 
Maisonobe's method for the s-parametric annihilator Ann^ n [ s ](/ S ) C -D„[s], where s = 
(si, . . . , s p ). The corresponding procedure in dmod.lib is called annfsBMI. It computes 
both Ann£ n [ s ] / s C D n [s] and the Bernstein-Sato ideal in K[s], which is defined as 

B(f) = (Ann Dn[si ,..., Sp] (A sl • • • i7) + </i--- U)) n K[ Sl , ...,s p ]. 

In contrary to the case / = f\, in general the ideal 13(f) need not be principal. How- 
ever, it is an open question to give a criterion for the principality of 13(f). Armed with 
such a criterion, one can apply a generalization of the method of Principal Intersection 
5.11 to multivariate subalgebras [2] and thus replace expensive elimination above by the 
computation of a minimal polynomial. Otherwise we still can apply the Principal Inter- 
section, which, however, will deliver only one polynomial to us. As in the case / = f± it is 
an open question, which strategy and which orderings should one use in the computation 
of the annihilator and of the Bernstein-Sato ideal in order to achieve better performance. 

We reported in [20] on several challenges, which have been solved with the help of our 
implementation. Namely, the products (x 3 + y 2 )(x 2 + y 3 ) and (x 2 + y 2 + y 3 )(x 3 + y 2 ) give 
rise to principal Bernstein-Sato ideals. 

Example 9.1. Let us consider the following example from [3], which is quite challenging 
to compute indeed. 

LIB "dmod.lib"; ring r = 0, (x,y , z) ,dp; 
ideal F = z, x~5 + y~5 + x~2*y~3*z; 
def A = annfsBMI (F) ; setring A; 
LD; // prints the annihilator in D[sl,s2] 
BS; // prints the Bernstein-Sato ideal 

We do not show the output here because of its size. But from the output one can see, 
the Bernstein-Sato ideal is of dimension 1 in /^[si,,^] since its Grobner basis consists 
of three elements {15625sis^ + 17 l.o.t., 3125s?4 + 24 lo.t., 625sfs^ + 42 l.o.t.}. Notably, 
every generator factorizes into linear factors and each factor involves either s\ or s%, which 
happens quite seldom in general. 

In general, quite a little is known about Bernstein-Sato ideals. Their dimensions, prin- 
cipality, factorization of generators and primary decomposition constitute open problems 
from the theoretical side. 
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10 Bernstein-Sato polynomial for a variety 

Now we proceed to the construction of the Bernstein-Sato polynomial of an affine algebraic 
variety. We refer to [7] for the details of the complete construction for arbitrary varieties 
and to [1] for the details about annihilator-driven algorithms. 

Given two positive integers n and r, for the rest of this section we fix the indices i, j, k, I 
ranging between 1 and r and an index m ranging between 1 and n.Let / = (fi, ■ ■ ■ , f r ) 
be an r-tuple in K[a;] r . Here, s = (si, . . . ,s r ), \ = , -_ , and f s = f^ 1 --- f£ r - Let us 
denote by K{S) the universal enveloping algebra U(#l r ), generated by the set of variables 
S = (sij), i,j = 1, . . . ,r, with su = s^ subject to relations [s„, s kl ] = 5 jk s a - S u s kj . We 
denote by D n (S) := D n cg) K K(S'), which is a G- algebra of Lie type by e. g. [22]. Then 
the free K[x, s, -^-module of rank one generated by the formal symbol f s has a natural 
structure of a left D n (S')-module: 

Sij • (G(s) ■ f s ) = s t ■ G(s + e j -e t )^-f s e K[x, s, h ■ f s , 

Ji J 

where G(s) G K[x, s, 4] and €j stands for the j'-th standard basis vector. 

Theorem 10.1 (Budur et al. [7]). For every r-tuple f = (fi,---,f r ) G K[x] r there 
exists a non-zero polynomial in one variable b(a) G K[cr] and r differential operators 
P x (S), ... , P r (S) G D n (S) such that 

J2 P ^ S )h ■ f S = b ^ + ' ' • + Sr) • / s G K[x, s, j] ■ f s . (2) 

fc=i ^ 

The Bernstein- Sato polynomial bf(o~) of / = (j\, . . . , f r ) is defined to be the monic 
polynomial of lowest degree in the variable a satisfying the equation (2). It turns out, that 
every root of the Bernstein-Sato polynomial is rational, as in the case of a hypersurface. 
Let / be the ideal generated by /i, . . . , f r and Z the affine algebraic variety associated 
with / in K n . Then it can be verified, that bf(o~) is independent of the choice of a system 
of generators of J, and moreover that bz{&) = bf(o~ — codimZ + 1) depends only on Z. 

Now, let us denote by AnnD n {s)(f s ) the left ideal of all elements P(S) G D n (S) such 
that P(S) • f s = 0. We call this ideal the annihilator of f s in D n (S). From the definition 
of the Bernstein-Sato polynomial it becomes clear that 

(Ann Dn(s) (r) + (/^ . . . , f r )) n K[ Sl + ■ ■ ■ + s r ] = (b f ( Sl + ... + s r )). 

Since the final intersection can be computed with Principal Intersection 5.11, the above 
formula provides an algorithm for computing the Bernstein-Sato polynomial of affine 
algebraic varieties, once we know a Grobner basis of the annihilator of f s in D n {S). 

Theorem 10.2. Let f = (/i, . . . , f r ) be an r-tuple in K.[x] r and D n (dt, S) the "K-algebra 
generated by D n , dt and S with the non- commutative relations of D n (S), described above 
and additional relations [s^, dt k ] = 5j k dti (dti commute mutually with the subalgebra D n ). 
Then the annihilator of f s in D n (S) can be expressed as follows: 



D n (dt, S) ( Sij + dtifj , d m + J2 ^dt k 

fe=l 
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1 < i, j '< r 

1 < m < n 



nD n (S). 



Note, that this result and its proof [1] can be presented as natural generalization of 
the algorithm for computing Anno„[ s ](/ s ) with the method of Briangon-Maisonobe (cf. 
Section 3). 

As Budur et al. point out [7, p. 794], the Bernstein-Sato polynomial for varieties 
coincides, up to shift of variables, with the 6-function in [34, p. 194], if the weight vector 
is chosen appropriately, see also [37] . Algorithms for computing the 6-function have been 
already discussed in Section 5, so the procedure bf ct Ideal can be immediately applied to 
this situation. Hence, like for the case of a hypersurface, we have two essentially different 
ways to compute Bernstein-Sato polynomials for varieties. The comparison of these two 
methods is the subject of further research. 

In the new SINGULAR library dmodvar .lib 1 , we present the implementations of the 
following algorithms 

Sannf sVar, which computes Aim.£) n (s){f s ) according to the Theorem 10.2, 

bf ctVarln, which computes bf(s± + . . . + s r ) using initial ideal approach, 

bf ctVarAnn, which computes bf(s 1 + . . . + s r ) using annihilator-driven approach. 

Example 10.3. Let TX = V(x% + y$, 2x x 1 + 3y%yi) C C 4 the tangent bundle of X = 
V(x 2 + y 3 ) C C 2 . Then the Bernstein-Sato polynomial of TX can be computed as follows: 

LIB "dmodvar. lib"; 

ring R = 0,(xO,xl,yO,yl),Dp; 

ideal F = xO~2+yCT3, 2*x0*xl+3*y0~2*yl ; 

bf ctVarAnn(F) ; // annihilator-driven approach 

// alternatiely , one can run 

bf ctVarln(F) ; // approach via initial ideal 

In both cases we obtain the polynomial 

b TX (a) = (a + 1) V + 1/3) V + 2/3) V + l/2)(cr + 5/6)(cr + 7/6). 

The annihilator ideal can be computed via executing, in addition to the first 3 lines of 
the above code, the following code: 

def S = Sannf sVar(F) ; // returns a ring 

setring S; //in this ring, ideal LD is the annihilator 

option(redSB) ; LD = groebner (LD) ; // reduced GB of LD 

There are 15 generators in the Grobner basis of AnnF s : 

3y 2 <9:ri - 2x dy 1 , 3y%dx + 6y yidx 1 - 2x dy - 2xidyi, xoyodx^yo - xoyodxody! + 
xiy dxidy 1 -2x yidxidyi, 3y y 1 dxl-x dx 1 dyo+x dx dy 1 -x 1 dxidy 1 , Sxoyoyidxodx^y^ 
6xoyfdxldyi — x^dxidy^ + x^dx^dy^dyi — 2x xidxidyodyi + xoXidxodyf — x\dx\dy\ + 
3xiy 9xi+3x yi9xf-3y 2/i9xi9yi, 6x yfdxjdyody 1 +3xoyoyidxldyf-3x 1 y y 1 dx dx 1 dyf+ 
6xoyfdxodxidyf—xldxidyQ+xldxodyldyi—2xoXidxidyQdyi+xoXidxodyodyf—x1dxidyodyf+ 
3x 1 y dx 2 1 dyo+3x yidx 2 1 dyo+9x y 1 dxodx 1 dy 1 -6y y 1 dx 1 dyodyi+3y y 1 dx dyf+6yfdx 1 dyl+ 
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3x\dx\ + 3yidxidyi, §xoy\dx\dyi — x^dx\dy^ + 2xldxodxidyodyi — 2xoXidx^dyodyi — 
xldxldyl+2xoXidxodxidyf—xldxldyf—3xoyodxodxi+3xiyodxl+3xoyidxl—2xodxidyodyi + 
x dx dyf-3x 1 dxidyf + 6yodxj, s 22 -x 1 dx 1 -y 1 dy 1 , 6s 21 -3x dx 1 -2y dy 1 , 6s u - 
3x dx +3x 1 dx 1 -2y dy +iy 1 dy 1 , s 12 x +3y yjdx 1 -x x 1 dx +x 2 1 dx 1 -x y 1 dyo, si 2 y dyi- 
x 1 yodxidyo + 2xiyidxidy 1 -yoy 1 dy dy 1 + 2yldyj-yody +4y 1 dy 1 , 3s 12 y 2 + 6x 1 y y 1 dx 1 - 
3y 2 y 1 dy +6y y 2 dy 1 -2x x 1 dy , s 12 x 1 dyf-3s 12 y dx 1 +3x 1 yoy 1 dx dx 1 dy 1 +6x 1 yldxjdy 1 - 
3y yfdx 1 dyody 1 +3yoyfdxodyf+6yfdx 1 dyf-x x 1 dx 1 dyo+xox 1 dxodyody 1 -2xldx 1 dy dy 1 - 
x 1 y 1 dy dyl+3x 1 yodx dx 1 +3x 1 y 1 dxl-6y y 1 dx 1 dyo+9yoyidx dy 1 +18yldx 1 dy 1 -2x 1 dy dy 1 + 
3y dx , 3s 12 y y 1 dx l -s 12 x l dyi-3x 1 y y 1 dx dx 1 -3y yfdx dy 1 +xjdx 1 dy +x 1 y 1 dy dy 1 - 
3y yidx + xidy . 

This ideal belongs to the K-algebra D^S) in 12 variables, as defined in the beginning of 
this section. By executing 

GKdim(LD); 

we obtain, that the Gel'fand-Kirillov dimension of D±(S)/ AnnF s is 6, the half of the 
Gel'fand-Kirillov dimension of D±(S). However, D±(S)/ AnnF s is not a generalized holo- 
nomic /^(S^-module, since the annihilator of this module contains a central element 
S12S21 — S11S22 — S11 and hence is not zero. 
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